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We prove that the predicted charge transfer state in symmetric bilayers of two dimensional electron 
gases is always unstable at zero bias voltage, due to interlayer correlation and/or tunneling. This 
is most easily seen by resorting to a pseudospin formalism and considering coherent states obtained 
from the charge transfer state through rotations of the pseudospins. Evidently, the charge transfer 
state is stabilized by a sufficiently strong gate voltage, as found in recent experiments. We show that 
a simple model, in which the layers are strictly two dimensional, is able to account quantitatively 
for such experimental findings, when correlation is properly included. 
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Two dimensional electron gas (2DEG) systems con- 
fined in semiconductor space charge layers (e.g., Si inver- 
sion layers, GaAs heterojunctions, and quantum wells) 
have provided for the last twenty years an ideal labora- 
tory for studying various electron-electron interaction ef- 
fects under almost ideal 2D jellium conditions. The zero 
temperature phase diagram has attracted a lot of atten- 
tion since Wigner pointed out that at low density elec- 
trons would crystallize to minimize the potential energy. 
A stable spin polarized phase at intermediate density was 
predicted by Hartree-Fock (HF) calculations, as a con- 
sequence of competition between kinetic and exchange 
energy. The most recent Diffusion Monte Carlo (DMC) 
simulations by Rapisarda and Senatore (RS) |l| confirm 
this picture and show that the 2DEG is paramagnetic 
up to r s — 20, then undergoes a ferromagnetic transition 
and eventually crystallizes at r s = 34. As usual, the cou- 
pling strength r s is defined in terms of the areal density 
by U2D = l/Tr^as) 2 , o,b being the Bohr radius. 

Bilayered systems, in addition to enhanced correlations 
effects with respect to single layers, have an additional 
degree of freedom — the interlayer distance d — and ex- 
hibit a more complex phase diagram. Paralleling a rich 
discussion on the high magnetic field Fractional Quan- 
tum Hall regime and references therein], some novel 
interaction-driven phases have recently been proposed for 
the zero field case [|H|-|j|- 

Based on the well-known fact that the compressibil- 
ity of the low density 2DEG is negative, Ruden and 
Wu (RW) [|| argued that exchange and correlation could 
overcome the kinetic and Hartree energy costs stabiliz- 
ing a charge transfer state (CTS), in which one layer 
contains all the electrons and the other is empty. Their 
HF computation, restricted to the ideal zero-tunneling 
symmetric case, predicts stability of the CTS for 
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tainable electron densities for reasonable layer separa- 
tions d (e.g., with GaAs parameters, r s = 4 or ri2D = 
2 ■ 10 10 cm~ 2 for d < 90 A, according to Ref. @). In 
fact higher values of the coupling (r s w 20) have been 
achieved [Io[ working with holes, rather than electrons. 
Note that in the bilayer we relate the density /coupling 
parameter r s to the total areal density. 

Evidence for a CTS in GaAs/AlAs double quantum 
well structures in the presence of bias voltage was re- 
cently reported |pd|-^3|] . The experimental setup allowed 
to measure independently the charge contained in the 
two layers as a function of bias voltage, via Shubnikov- 
de Haas oscillations. For low voltages the transferred 
charge was found to depend linearly on V e , in quantita- 
tive agreement with the predictions of a simple equivalent 
circuit model, but for high enough voltages V g ^ IV the 
charge of one layer dropped abruptly from around 20% 
of the total charge to 0. 

Here we shall demonstrate that, for two dimensional 
electron layers and at zero bias voltage, the CTS is un- 
stable at any density and layer separation d. A conve- 
nient model Hamiltonian for the coupled electron layers 
is obtained introducing a pseudospin variable n, which 
labels the planes with the convention that t- = 1/2 
(—1/2) if electron i is in the upper (lower) plane. For 
fixed interlayer distance and background densities and in 
the absence of tunneling, the Hamiltonian takes the form 
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which is well into the range of the experimentally at- 



where V e is the energy difference between the two planes, 
given by 2ire 2 (pi — p2)d in the absence of external fields, 
pi and P2 being the background densities of the two 
planes. The prime on the second sum indicates that 
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the term with k = should be omitted when the sum 
is rewritten in reciprocal space. It is clear that with 
N electrons there are two CTS's, corresponding to two 
eigenstates of the total pseudospin component along z, 

T(z) = J2i A Z) = ±N/2. In other words CTS is the same 
as (pseudospin) ferromagnetism, with the magnetization 
oriented along the z axis. We shall show that without 
bias voltage there can be no such ferromagnetic phase 
transition. If any, in the presence of tunneling a transi- 
tion takes place to a (coherent) correlated state, which 
has the same symmetry as the charge transfer state but 
does not displace any charge. The charge transfer phase 
transitions can take place in the presence of significant 
sample asymmetry and/or gate voltage. 

It is interesting to consider the ideal case d — V e = 0, 
where both real spin a and pseudospin r are conserved 
and the system is actually a 4-component monolayer 
2DEG. At large density (small r s ) the ground state 
is completely paramagnetic, with the four components 
equally populated to minimize the Fermi energy; in the 
opposite low density regime the ground state is known to 
be a Wigner crystal and essentially insensitive to spin po- 
larization. In the intermediate regime the 4-component 
2DEG can be expected to mimic the above-mentioned 
behavior of the 2DEG (and 3DEG), with the appear- 
ance of a spin-polarized ground state. In HF, the 4- 
component 2DEG is unstable towards the 2-component 
phase at r s — 37r/16(-\/2 — 1) ~ 1.42 and then towards 
the 1-component phase at r s = 3n/8V2{V2- 1) ~ 2.01, 
the first being obviously identical with the d — > limit 
of the RW results and the latter with the usual spin po- 
larization transition of the 2DEG in HF. 

To assess the actual presence of these phase transitions 
we performed Slater-Jastrow Variational and Fixed Node 
Diffusion Monte Carlo simulations for the 4-component 
2DEG, and compared them with the RS results. The 
method is completely analogous to the one employed by 
RS, with the only difference that our DMC code moves 
all particles at every step, therefore requiring smaller 
time steps which, in turn, makes time step extrapola- 
tions unnecessary. To check the code we duplicated one 
of the RS data points, finding excellent agreement within 
statistical noise. Our results are reported in Table Q, 
and compared in Figure [l] with the data by RS on the 
2 and 1-component 2DEG. At variance with the sim- 
ple HF prediction, the figure clearly shows that the 4- 
component phase crystallizes at r s ~ 42, being always 
stable at smaller couplings with respect to the other two 
fluid phases that we have considered. One could argue 
that backflow corrections might spoil the validity of our 
results, but we do not think this to be the case since 
backflow corrections are known to lower the paramag- 
netic energies more than the ferromagnetic ones. 

As we have already mentioned, a total charge trans- 
fer state is described in our formalism as an eigen- 
state of each with eigenvalue +1/2. We remark 

(z) 

that t- commutes with the Hamiltonian even at fi- 



FIG. 1. Ground state energy as a function of r 3 for the nor- 
mal unpolarized (E2), polarized (Ei) and 4-component (E4) 
2DEG and for the triangular Wigner crystal (E so i). Except 
for Ei, all data are from RS. On the purpose of clarity we 
plotted r 3 J 2 (E{r s ) - Ci/r s ), with ci = -2.2122. 
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nite d and V e (but zero tunneling) and we label with 
\^ ch.tr) the lowest eigenstate within the charge trans- 
fer subspace. In the symmetric V e — case it 
can be easily seen that the coherent state \^ co h) — 

exp |«2 J2j T j^} \^ ch.tr.) has an average energy lower 
than \^ ch.tr.)- In fact, (i) l/x > l/y/d 2 + x 2 for 
any x; and (ii) ch.tr.^T^ t^ z) \^> ch . tr ) = 1, whereas 

coh.l^T^ Tj \^ coh.) =0- The variational principle im- 
plies that \^ ch.tr.) cannot be the ground state of the sys- 
tem, even if \^coh.) is n °t an eigenstate of H. 

We stress that both \^ c h.tr.) and \^ C ah.) have the same 
symmetry properties under exchange of the particle's 
spatial coordinates, as this is not affected by pseudospin 
rotation. In particular the two states are both fully pseu- 
dospin polarized, though in different directions. The av- 
erage charge distribution in \^ co h.) corresponds to half 
electrons in one layer and half in the other, and therefore 
it does not loose in Hartree energy with respect to the 
normal state, still having the exchange energy gain deriv- 
ing from the full antisymmetry of the spatial part of the 
wavefunction. The difference with the above-mentioned 
HF predictions is due to the fact that no state with the 
same symmetry as \&coh.) was considered by RW. 

No such clearcut statement can be made in the asym- 
metric V e 7^ case. In general, it can be seen that the 
term V e J2i favours finite values of , i.e. (partial) 
pseudospin polarization. A quantitative estimate would 
rely on the detailed pseudospin dependence of the corre- 
lation energy for the bilayer system, which — to the best 
of our knowledge — has not yet been determined. Com- 
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mon wisdom saying that the phase transition in jellium 
models occurs directly from unpolarized to fully polar- 
ized implies the presence of a maximum in the internal 
energy at intermediate polarizations. In turn, this im- 
plies that if the value of V e is sufficient to raise the po- 
larization up to this point, which may still be far below 
saturation, it will spontaneously saturate. This behav- 
ior would look like a 2-component to 1-component phase 
transition driven by an increasing V e , and can poss ibly 
explain the experimental results reported in Refs. |Tl| . |l2| . 

It must be noticed that the state \^ C oh.) is a broken 
symmetry state, since all states of the form co h.{4>)) = 
exp |«| Sj( T j !/) cos0 + t^ x) sin^j \^ ch . tr ) are degener- 
ate. The order parameter is the total pseudospin T = 
J2i T i> corresponding to interlayer phase coherence, and 
behaves as an easy-plane ferromagnet, the z component 
being frozen by the interplay of V e and the Hartree en- 
ergy. Most of the considerations presented in Ref. Q for 
the high B situation can be directly extended to this zero 
field case. The role of a tunneling Hamiltonian, which can 
be conveniently written as 

H t = tJ2rt\ (3) 

i 

with the tunneling matrix element t > , is to break the 
(j) symmetry and stabilize the coherent state \^ co h.) = 
\V coh .((l) = 0)), in which (T< z )) = and = -N/2. 
In the large t case all electrons trivially lie in the sym- 
metric state, which has exactly the same form as \^ C oh)- 
For small enough layer separation d and tunneling ampli- 
tude t a perturbative estimate of the critical amplitude 
t* at which the fully coherent state \^ C oh) becomes sta- 
ble is easily obtained from the energies of Table || and of 
RS, on the ground that the d dependence is negligible, 
as it vanishes in first order, contrary to the tunneling 
term. Thus one gets f*/2 = Ei(r s ) — E 2 (r s ), which at 
r s = 2 corresponds to t* — 0.8 meV = 9K with GaAs 
parameters. 

The experimental setup of Ref. |Tl|-[T3]] includes a 
metallic gate at a large distance D from the two 
2DEG's and some charged dopants, whose location and 
amount are unspecified and irrelevant for the present pur- 
poses. The 2DEG's are realized with 150(180)A thick 
GaAs quantum wells with a midpoint separation d of 
220(194)A for sample A(B), and the tunneling energy 
t = 0.005K (5K). From the work of Zheng and Mac- 
Donald |H| and from numerical simulations B it is clear 
that interlayer correlations are relevant only if d becomes 
smaller than the average interparticle distance r s as- As 
in the present situation r s aB ~ 150 A and d ~ 200 A cor- 
relation effects can be safely neglected. In the following 
we shall further neglect finite thickness and tunneling ef- 
fects. 

We propose a simple model which correctly includes 
intralayer correlations and interlayer mean field interac- 
tions and is closely related to the one used by Eisenstein 
et al pa] to discuss a similar experimental setup. Let us 



FIG. 2. Layer densities for the bilayer A specified in the 
text, as a function of the bias voltage. Triangles and diamonds 
give the experimental results of Ref. [Q. The full curve re- 
sults from the strictly two dimensional model discussed in 
the text, while the dashed curve is obtained neglecting both 
correlation and exchange. The inset reports the pseudospin 
polarization £ as a function of bias potential V g . 
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consider the gate potential Vg corresponding to equal 
density in the two 'active' layers, i.e. to zero elec- 
tric field between the layers. [The value of Vg ^ depends 
on details of the sample design, as location and amount 
of dopants, geometry, etc.]. With reference to this state 
and under the assumption that no net charge flows in the 
system with varying V g , the electrochemical equilibrium 
condition between the two 2DEG's is written as 

H(rn) - M(n (0) ) = fi{n 2 ) - M (n (0) ) + Ane 2 d(n 2 - n ) (4) 

and between the upper layer and the gate as 

- e(V g - V g {0} ) = /z(m) - ^(n (0) ) + 47re 2 L>(ni + n 2 - 2n (0) ) 

(5) 

We now proceed to solve numerically these two coupled 
equations for the unknowns n\ and n 2 at various values 
of V g . While the sample parameters d, and V^ -* are 
given in Ref. |]l2|| , D is not, and was determined by fitting 
the measured variation of the average density (rii +n 2 )/2 

with Vg. 

The chemical potential fjt(n) = d(nE(n))/dn for the 
idealized single layer 2DEG is derived from the precise 
equation of state obtained by RS with DMC. The re- 
sults shown in Figs. ||, || are in very good agreement 
with the experimental data from Ref. jl2] ], with signif- 
icant disagreement only in the tunneling-dominated re- 
gion around rii = n 2 . In fact quantitative agreement (see 
Fig. ||) can be obtained even in this region computing the 
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FIG. 3. Layer densities for the bilayer B specified in the 
text, as a function of the bias voltage. Triangles and diamonds 
give the experimental results of Ref. [Q. The full curve re- 
sults from the strictly two dimensional model discussed in 
the text, while the dashed curve is obtained by accounting 
for tunneling as explained in the text. 
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subband densities from the layer densities with no tun- 
neling hi, fi2 by exact diagonalization of a noninteracting 
electron hamiltonian with finite tunneling, 
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M°(«2) j 4,2 a fe,2 + t («L a fe,2 + CC.) 
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We remark that neglecting intralayer interactions would 
lead to the usual approximation fjP(n) — C q n 7 where the 
so-called quantum capacitance C q is given by irfi 2 jm. By 
making this assumption for both layers the above equa- 
tions are linear and can be solved analytically, obtaining 
results in fair agreement with experiments but completely 
missing the nontrivial charge transfer effect signaled by 
the peak in the upper curve of Fig j^. 

In conclusion, we presented a pseudospin formalism to 
describe bilayer 2DEG systems and used it to rigorously 
prove the instability of the charge transfer state and to 
argue for the stabilization of coherent states at small fi- 
nite tunneling. We also reported results from DMC sim- 
ulations, which reveal the stability of the pseudospin un- 
polarized state at d = V e = 0. Finally we have shown 
that the presently available experimental data, being in 
a regime where interlayer correlations are negligible, can 
be quantitatively accounted for by a simple strictly two 
dimensional model. Clearly the more challenging and in- 
teresting situations, in which interlayer correlations are 
discernible and determine the details of the charge trans- 
fer, is still awaiting for both experimental and theoretical 



investigations. 

A preliminary account of the above study was pre- 
sented at the INFM-FORUM workshop on 2DEG, held 
in Pisa, FORUM, June 96, and at the Conference The 
Electron Quantum Liquid in Systems of Reduced Dimen- 
sions held in Trieste, ICTP, July 96. We acknowledge 
stimulating discussions with B. I. Halperin and with A. 
H. MacDonald. 

Note added — While completing the present manuscript 
we became aware of a recent preprint on the same topic 
by Das Sarma et al |l6| , which is somewhat complemen- 
tary to this work. 
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TABLE I. Fixed-node DMC total energies of the 4-component 2DEG for 52 particles and in the 
bulk limit, in Ryberg per particle. Variational results used for the size-extrapolation are also given. 
These data are accurately reproduced by the same fitting formula used by RS, with parameters 
ao = -0.88115, 01 = 4.439, a 2 = 0.14063, a 3 = 1.9034. 





r s = 2 


r s = 10 


r s = 20 


r s = 30 


r s = 50 


N = 36 


-0.5882(4) 


-0.17102(2) 


-0.092275(6) 


-0.063561(5) 


-0.039386(3) 


N = 52 


-0.5728(1) 


-0.17030(2) 


-0.092085(4) 


-0.063474(3) 


-0.039354(1) 


N = 84 


-0.58062(7) 


-0.17055(2) 


-0.092123(5) 


-0.063481(4) 


-0.039355(2) 


N = 100 


-0.5758(1) 


-0.17036(2) 


-0.092062(4) 


-0.063457(4) 


-0.039344(2) 


N = 180 


-0.57676(8) 


-0.17031(2) 


-0.092050(5) 


-0.063453(4) 


-0.039340(2) 


VMCvo 


-0.5750(1) 


-0.17019(3) 


-0.092002(4) 


-0.063426(3) 


-0.039330(1) 


DMC52 


-0.5838(4) 


-0.17232(4) 


-0.092891(9) 


-0.063975(8) 


-0.039639(3) 


DMC X 


-0.5860(5) 


-0.17221(6) 


-0.09281(1) 


-0.063927(9) 


-0.039615(4) 
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